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The 2-10 keV emission properties of PSR B1937+21 
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Abstract. We present the results of a BeppoSAX obser- 
vation of the fastest pulsar known: PSR B1937+21. The 
- 200 ks observation (78.5 (34) ks MEGS (LEGS) expo- 
sure times) allowed us to investigate with high statistical 
significance both the spectral properties and the pulse pro- 
file shape. The absorbed power law spectral model gave a 
photon index of 1.7 and ~ 2.3 x 10^^ cni-2. These 
values explain both a. the ROSAT non-detection and b. 
the deviant estimate of a photon index of ~ 0.8 obtained 
by ASGA. The pulse profile appears, for the first time, 
clearly double peaked with the main component much 
stronger than the other. The statistical significance is lOcr 
(main peak) and 5cr (secondary peak). The 1.6-10 keV 
pulsed fraction is consistent with 100%; only in the 1.6-4 
keV band there is a ~ 2cr indication for a DC component. 
The secondary peak is detected significantly only for en- 
ergies above 3-^4 keV. The unabsorbed (2-10 keV) flux 
is F2-10 = 3.7 X 10"^'^ erg cm"^ s~^, implying a lumi- 
nosity of Lx = 4.6 X 10^* 9 (d/3.6 kpc)^ erg s~* and an 
X-ray efficiency of 77 = 4 x 10~^ O, where G is the solid 
angle spanned by the emission beam. These results are in 
agreement with those obtained by ASCA. 



1. Introduction 

The various X-ray satellites, from Rosat to Chandra, 
observed about 10 millisecond pulsars (MSPs) (see e.g. 
Becker 2001), but good spectral and temporal informa- 
tions exist only for half of them. In fact MSP X-ray obser- 
vations are usually affected not only by low statistics but 
also by insufficient time accuracies to perform detailed pe- 
riodicity and timing analysis. X-ray observations demon- 
strate that the X-ray emission from MSPs is generally not 
of thermal origin, and those with the hardest spectra ap- 
pear to be objects with strong magnetic fields at the light- 
cyfinder Bl [Rl = cP/2it; see Saito et al. 1997, Kuiper et 
al. 1998, Kuiper et al. 2000, Becker 2001). Given the cor- 
relation between spin-down power and X-ray luminosity 
(see Verbunt et al. 1996, Becker and Triimper 1997 and 
Takahashi et al. 2001), another important quantity is the 
spin-down flux density at Earth {E/Att(P). 



PSR B1937-f-21, with a period P = 1.56 ms, was the 
first and still the fastest MSP known. In spite of its low 
surface magnetic field strength of Bs = 4.1 x 10^ G, its 
magnetic field at the light-cylinder is the highest of all 
known pulsars: Bl ^ I x 10^ G, very similar to that of 
the Crab pulsar. Its spin-down flux density is ~ 2 x 10""'^° 
erg cm"2 s~*, which is 3 times higher than for PSR 
J0218-F4232 but 3 times lower than for PSR B1821-24. 

Here we report on the results obtained by a BeppoSAH 
observation in the energy range 1.6-10 keV. 

2. Observation and spatial analysis 

BeppoSAX observed PSR B1937+21 on May 1"* 2001. The 
~ 200 ks observation time resulted in a total effective 
on-source time of ^ 78.5 ks (~ 34 ks) with the MEGS 
(LEGS). To extract the source photons from LEGS (0.1- 
10 keV) and MEGS (1.6-10 keV), we used: 

— a standard spatial analysis with optimized extraction 
radius and energy range; 

— a Maximum Likelihood (ML) approach to extract the 
number of counts assigned to the pulsar taking into 
account the presence of other sources and the back- 
ground simultaneously. 

In the first case we estimated the local background (in a 
circular region in the field of view 10' away from the pul- 
sar) and compared the results with the standard one from 
archival blank sky usually used in the spectral fitting. We 
found that the local background is 50% and 15% higher 
for LEGS and MEGS respectively. The extraction radii 
which optimized the signal to noise ratio were 3' (LEGS) 
and 2' (MEGS) resulting in a collection of 106 (0.5-8 keV) 
and 369 (1.6-10 keV) photons. In spite of the lower expo- 
sure time and sensitivity of LEGS we were able to detect 
a significant signal in the 0.5-8 keV range for this instru- 
ment. Inclusion of these softer photons in the spectral fit 
allowed us to better constrain the hydrogen column den- 
sity toward the pulsar (iVn)- Using the ML approach we 
could extract the source photons modeling the contribu- 
tion to the total counts from additional field sources and 
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the (assumed flat) background. This method was success- 
fully used in several other cases of weak source (see e.g. 
Kuiper et al. 1998, Mineo et al. 2000). 

3. Temporal analysis 

The arrival times of the extracted photons were converted 
to the Solar System Barycentric (SSB) Frame using the 
BARYCONV]^ code and then searched for periodicities by 
folding them with trial frequencies around the radio one 
(see updated radio ephemeris in Table |l|). The value 
plotted versus trial frequency is shown in Fig. ^ The dis- 
tribution had a clear maximum, but we found the corre- 
sponding pulse frequency deviating from the radio value 
by ~ -7.6 X 10^^ Hz (/;,; — fx)- This is about ten times 
the statistical error of ~ 8 x 10^^ Hz. To investigate this 
discrepancy we first verified the correctness of the SSB 
conversion code by also using independent software. We 
verified that all the routines involved in the processing 
were working properly. We then verified the OBT (On- 
Board Time) to UTC time conversion routine. To this aim 
we: 

1. analyzed separately two sections of the observation; 
the results confirmed those of the full observation; 

2. used different polynomial order for the OBT-UTC fit 
and various high deviation points (de-)selection crite- 
ria; in all cases the results were the same; 

3. verified the (X-ray) pulse frequency of the Crab dur- 
ing the observation performed on April 2000 and on 
September 2001; these observations arc shorter than 
the PSR B1937+21 one (in particular the 2001 ob- 
servation), nevertheless in the first case wc found a 
discrepancy which is outside the statistical errors. 

We concluded that the difference had likely to be at- 
tributed to incorrect photon time markers. We investi- 
gated the possible reasons of that with the BeppoSAX 
mission director. The answer was that, excluding an un- 
likely (and unobserved) UTC clock drift (which comes out 
of the Malindi ground station GPS clock), the only likely 
explanation could be an on-board electronics problem. In 
particular it could be an internal drift of the clock which 
tags the "base band" packets. Performance degradation 
with time of the clock is not unexpected and usually not 
considered critical for detectors not specifically built for 
high precision timing. Also, as we verified with the Crab, 
the problem is less severe (or not visible) for pulsar with 
relatively long periods or for relatively short observations. 
Even though other (we believe more unlikely) explanations 
could be found, we decided to correct the (possible) clock 
drift by estimating the required additional slope to add 
in the OBT-UTC conversion. The result was a constant 
slope of ~ 1 X 10~* s s~^. We verified this value (actually 
a value slightly lower) could also correct the 2000 Crab 



Table 1. Ephemeris of PSR B1937+21 (from observations 



performed at Effelsbcrg). 


Parameter 


Value 


Right Ascension (J2000) 


19'^ 39'" 38=560 


Declination (J2000) 


21° 34' 59'.' 14 


Epoch validity start/end (MJD) 


51639 - 52047 


Frequency 


641.928246349481 Hz 


Frequency derivative 


-4.331 X IQ-" Hz s-^ 


Epoch of the period (MJD) 


51843.0 


RA proper motion 


—0.15 mas yr~^ 


Dec proper motion 


—0.47 mas yr"^ 



140 - 
120 - 
100 - 




-IxlO ^ -5x10 ^ 5x10 ^ IxlO 



Frequency - F„ (Hz) 

Fig. 1. The significance as function of trial folding 
frequency. The is aligned with the maximum. 

observation ('^ 70 ks elapsed time) while it had no effect 
on the 2001 observation 24 ks elapsed time). 

Counting statistics allowed us to perform temporal 
analysis only for the MECS data (actually a search on the 
LECS photons didn't reveal any significant modulation). 
As already mentioned, a periodicity search analysis was 
performed on the 369 photons (1.6-10 keV) around the 
nominal radio period, after converting the arrival times to 
the SSB. Figure ^ shows the folded light curve together 
with the radio one (from the EPN archive^ - Kramer et 
al. 1998). The primary peak (PI) is detected at lOcr level 
and, for the first time, wc also detect the secondary peak 
(P2 - significance ^ 5cr) and estimate a relative phase 
separation of 0.48 ± 0.04. Data collected by ASCA al- 
lowed Takahashi et al. (2001) to claim an alignment of 
the main X-ray peak with the radio intcrpulsc. Wc then 
arbitrarily aligned the two profiles (the BSAX on-board 
clock absolute calibration docs not allow to do it for a 



^ see http://www.sdc.asi.it/software/saxdas/baryconv.html 



^ see http://www.mpifr-bonn.mpg.de/pulsar/dat/ 
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Table 2. Power law fit parameters using photons ex- 
tracted with the standard and the Maximum Likehhood 
methods. 



Method 


AT" 


a 




F2 


b 

-10 








(10^2 cm-2) 






STD 


0.85 


1-OO-0.06 


9 97+0.58 
^■^ ' -0.47 


3.7 


±0.4 


ML 


1 19+0.10 
'-■'-'^-0.10 


1 71+0.06 
' -^-0.08 


2 ^o+°■^^ 

^■l"-0.42 


3.7 


±0.4 



Note: all quoted uncertainties are la confidence. It is 
NphiE) = KxE-°'. 

Normalization at 1 keV (10^ ph s~^ cm~^). 
Unabsorbed flux in 2-10 keV (10"^^ erg cm~^ s"^). 



MSP). The vertical dashed lines in the figure are exactly 
0.5 apart. The radio peaks separation is instead 0.4794(8) 
so it would be tempting to consider the ASCA result in- 
correct and consider PI be coincident with the main radio 
peak. The X-ray peak widths are compatible with the in- 
strument time resolution limit (~ 100 /xs). Figure ^ shows 
the X-ray pulse profile in 3 energy ranges. In an attempt 
to quantify the evidence for a pulsar DC component, we 
applied the bootstrap method proposed by Swanepoel et 
al. (1996), which allows us to estimate the DC level (sky 
background and possibly a source component) in the pulse 
profile. We can then determine the DC and pulsed frac- 
tions of the pulsar by taking the background level from 
the spatial analysis into account. However, in the available 
form, the bootstrap method is able to find only one un- 
pulsed interval and cannot account for systematic errors. 
One should therefore realize that the quoted errors on the 
calculated parameters are only statistical. Note in Fig. || 
how the DC level is compatible with the background level. 
We find that over the total energy range 1.6-10 keV the 
source is consistent with being 100% pulsed; the 2a lower 
limits for the pulsed fraction in the energy ranges 1.6-10, 
1.6-4 and 4-10 keV are 88%, 66% and 95%, respectively. 
Only for the lower 1.6-4 keV energy interval we find an 
indication for a DC component of (18.7 ± 7.7)%. An in- 
crease of P2 compared to PI with energy is clearly visible. 



4. Spectral analysis 

The spectral analysis was performed using photons ex- 
tracted in two different ways: a. those in a (S/N optimized) 
circular region centered on the pulsar position (standard 
method); b. those from the ML analysis, which consider 
simultaneously the source, the background and a number 
of field sources. In both cases the best fitting model was 
an absorbed power law. A Black Body model did not fit 
the data satisfactorily. The two methods gave consistent 
results and are summarized in Tab. The absorbed fiux 
(2-10 keV) is 3.1 x 10"" erg cm'^ 8"^ 




1.0 

Pulse phase 

Fig. 2. The 1.6-10 keV X-ray profile and the radio one 
arbitrarily aligned using the main X-ray peak. The two 
vertical dotted lines are 0.5 apart. The X-ray phase sep- 
aration is 0.48 like for the radio but 0.5 shifted in phase 
(if the ASCA result is correct). In this case we can also 
confidently exclude that the secondary X-ray peak is as- 
sociated with the feature visible in the main radio peak 
(which instead would mimic the X-ray feature in PI). 



The derived (unabsorbed) luminosity is Lx — 4.6 x 
10'^^ O (c?/3.6 kpc)^ erg s~^ while the X-ray conversion 
efficiency is 77 = Lx/E = 4 x 10^^ 8, where 8 is the 
solid angle spanned by the emission beam. Note that the 
column density is 10 times greater than the radio DM de- 
rived one (by assuming 10 H atoms for each e~) and is 
only marginally consistent (within 2a) with the Galactic 
value 1.4 X lO'^^ cm~^. In addition the photon index is 
significantly softer than that found for PSR J0218+4232 
(0.94±0.22) and PSR B1821-24 (~ 1.1 using RXTE data, 
Kawai and Saito 1999). Also for the latter pulsar (placed 
in the globular cluster M28, at a distance of 5 kpc, 
galactic latitude — 5?6) its TVh — 3 x 10^^ cm^^, consis- 
tent both with the Galactic and DM derived one. Being 
the fine of sight to PSR B1937±21 tangent to the Galac- 
tic spiral arms {I ~ 57.51, b ~ —0.29, see Taylor and 
Cordes 1993) one can invoke a particularly low e^/iVn ra- 
tio in that direction due to the lack of ionizing sources. 
However, using the relation A^h = 1-79 x lO^^^y of Pre- 
dchl and Schmitt (1995), Verbunt at al. (1996) set an up- 
per limit of 2.2 x 10^^ cm~^ for the column density, i.e. 
in agreement with the DM value. More simply, the mea- 
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Fig. 3. The X-ray pulse profile in the full MECS X-ray 
band 1.6-10 keV and two sub-bands. The shaded areas 
show the estimated DC level (±lcr), while the dashed lines 
indicate the measured background level. In the 1.6-4 keV 
energy band the secondary peak is detected at <, 3a level 
and is compatible with the DC level. In the 4-10 keV band 
instead the detection level is 4(t (above background or DC 
level) . 



surcd high A^h could suggest for PSR B1937+21 a true 
distance greater than the 3.6 kpc derived from the Taylor 
and Cordcs (1993) model. If instead this value is correct 
then the continuing radio timing measurements should be 
able to detect its effect on the timing fit residuals very 
soon. We also note that the fit parameters are consistent 
with the non-detection in the ROSAT HRI energy range 
(0.1-2.4 keV). 

5. Conclusions 

We detected the double peak profile of the fastest pul- 
sar known. The 1.6-10 keV pulsed fraction is consistent 
with 100%, only for the band 1.6-4 keV there is an indica- 
tion for a DC component (18.7 ± 7.7 %). The secondary 
(X-ray) peak is significantly detected only above 3-4 keV 
and the ratio primary/secondary decreases with energy. It 
was not possible to perform an absolute phase comparison 
with the radio profile, but, purely morphologically, a phase 
separation comparison suggests alignment of the main X- 
ray peak with the main radio one, contrary to the ASCA 
finding by Takahashi et al (2001) which claim an align- 
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Fig. 4. The 1.6-10 keV spectrum and best fit model 
with the unabsorbed power law (a = 1.66). It is evi- 
dent the strong absorption below 3 keV. In the inset 
the smoothed MECS image around the PSR position is 
shown. 

ment with the radio interpulse. An accepted Rossi- XTE 
observation will clarify this point. 

We measure a soft power law spectral index {a ~ 1.7), 
different from those measured for PSR ,10218+4232 and 
PSR B1821-24. The Nu ^ 2.3 x 10^^ cm"^ marginally 
consistent with the Galactic value in the direction of the 
pulsar and gives a e~ /Nh ratio of 1/100, 10 times smaller 
than the canonical average value. This could be justified 
by the position of PSR B1937-I-21 in the Galaxy w.r.t. our 
line of sight but could also suggest a larger distance than 
the DM derived one. Radio timing should settle this issue 
soon. 
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